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Abstract 



f-H . It has recently been shown that identical, isotropic particles can form complex 

O ! crystals and quasicrystals. In order to understand the relation between the particle 

^ ' interaction and the structure, which it stabilizes, the phase behavior of a class of 

two-scale potentials is studied. In two dimensions, the phase diagram features many 
phases previously observed in experiment and simulation. The three-dimensional 
l/^ I system includes the sigma phase with 30 particles per unit cell, not grown in sim- 

r^ ■ ulations before, and an amorphous state, which we found impossible to crystallize 

in molecular dynamics. We suggest that the appearance of structural complexity in 
monodisperse systems is related to competing nearest neighbor distances and discuss 
(^ • implications of our result for the self-assembly of macromolecules. 

oo ■ 

p 

1 Introduction 

X 

^ [ When kept long enough at low temperatures, most systems develop long-range order. The 

easiest crystallization is expected for monodisperse systems, because only topological, but 
no additional chemical ordering is necessary. If the particles are isotropic, then a first 
reasoning suggests a preference to form simple, close-packed crystals, since this allows all 
of them to have the same, high first coordination numbers. A look at the periodic table 
reveals that indeed the ground states of most metals are bcc, fee, or hep [1]. Similar simple 
crystals are found in mesoscopic or macroscopic systems like for example globular proteins 
[2], monodispersed colloids [3], and bubble rafts on liquid surfaces [1]. 

In some of the systems above, quantum mechanics does not play an important role. 
Hence, it should in principle be possible to understand crystallization by using classical 
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pair interactions. For monodisperse systems there is only one type of interaction given 
by the shape of the potential function. Many common potentials are smooth and have a 
single minimum only. Simulations with these potentials usually lead to simple crystals. 
An example is the Lennard- Jones (LJ) potential. 

Nevertheless, the situation is not always that easy. Over the last years, more complex 
structures have been observed in experiments and simulations. A large part of work in this 
direction started after the discovery of quasicrystals in metallic alloys and was carried out 
with the aim to understand their formation. Another type of complex order are periodic 
crystals with large unit cells, so-called complex crystals. Whereas the lattice constants of 
a simple crystal are comparable to the range of the particle interactions, the unit cell of a 
complex crystal is stabilized indirectly, e.g. by geometric constraints. 

On the theoretical side, the obvious procedure to promote structural complexity is to 
use potentials with a more complicated radial dependence. It has been shown with double- 
minima potentials [5], oscillating potentials [6j, and a repulsive barrier [7j that the energy 
of an icosahedral phase can be lower than the energy of a class of trial structures including 
close-packed phases. Surprisingly, even in the LJ system, a quasicrystal is unstable only 
by a small energy difference [8]. 

These early works relied on general arguments and did not try to observe real crystals in 
simulations. Dzugutov was the first to actually grow a stable one-component dodecagonal 
quasicrystal from the melt |^, although it was later found to be only metastable [lO]. He 
used a LJ potential with an additional bump to disfavor the formation of simple crystals. 
Although no further work on other three-dimensional systems has been reported, there 
are many investigations focusing on two-dimensional systems, where computation and vi- 
sualization is easier. The first such study applied a variation of the Dzugutov potential 
and found a planar dodecagonal quasicrystal [TT]. Furthermore, it was pointed out that a 
decagonal quasicrystal appears with a square- well [12] and a ramp potential [T3] . 

Experiments show that nature is quite ingenious in her way to assemble identical par- 
ticles. First of all, a few complex ground states of elemental metals are known to exist. 
A notable example is a-manganese, which has cubic symmetry with 58 atoms per unit 
cell. Thermodynamically stable high-temperature phases are /3-boron with 105 atoms and 
/?-uranium with 30 atoms per unit cell [T]. The latter is isostructural to cr-CrFe and known 
as the sigma phase. Furthermore, commensurately modulated phases are common at high 
pressures p!^. 

Recently, different kinds of macromolecules have been observed to self-assemble into 
complex phases: (i) Under appropriate experimental conditions, tree-like molecules (den- 
drons) forming spherical micelles arrange to a dodecagonal quasicrystal [12] . (ii) T-shaped 
molecules can be designed in such a way that they organize into liquid crystalline honey- 
combs [16]. (iii) ABC-star polymers form cylindrical columns according to square-triangle 
quasicrystals, a two-dimensional version of the sigma phase ^7\ , and other Archimedean 
tilings. 

How can structural complexity be understood from bottom-up? It is instructive to 
study the ground state for small portions of the system - more or less spherical clusters - 
as a function of the particle number A^. The structure of small clusters (A^ < 200 in three 



dimensions, A^ < 30 in two dimensions) is often different from the bulk crystal [18]. The 
reason is the competition between local lowest energy configurations and the necessity of 
continuation in space. For example in a monodisperse LJ system, icosahedral coordination 
occurs in small clusters pTO], although the hep phase is the lowest energy bulk crystal ^20] . 
On the other hand, if small clusters already have simple structure, then the bulk phase 
will also be simple. A local order, which is incompatible with periodicity, is a necessary 
condition for structural complexity in the bulk. 

2 Phase diagram for the two-dimensional system 

There are two mechanisms to introduce structural complexity: either by destabilizing 
simple phases or by stabilizing a complex phase. An example for the destabilization mech- 
anism is the Dzugutov potential. We adopt the stabilization mechanism, because it has 
the advantage that the choice of the target structure can be controlled more directly. 

The particles are assumed to interact with an isotropic two-scale potential. A simple 
ansatz is the Lennard- Jones-Gauss (LJG) potential 

which has for small tq the shape of a shoulder, and otherwise represents a double-well 
with first minimum at r «i 1 with depth V{1) ~ 1 and second minimum around tq with 
V{ro) ~ e. The parameter a specifies the width of the second minimum. The shoulder 
form has been used extensively to observe liquid-state anomalies |121J and for understanding 
liquid-liquid transitions |22j. Surprisingly, almost nothing has been known until recently 
[231 [21] about the solid-state behavior of ([T]). 

In two dimensions, the local lowest energy configurations are regular polygons. As 
shown in Fig. [T](a), there are two possibilities: polygons with or without central atom. An 
n-gon without central atom is stabilized by having the second minimum at 

Vq = 2cos(7r/n). (2) 

The m-gon with central atom is favorable for 

rQ^ = 2sm{7r/m). (3) 

The values n = 3 and m = 6 lead to hexagonal local order (Hex), the value n = 4 to 
squares (Sqa), and n = 6 to honeycombs (Hon). Furthermore, pentagonal {n = 5, Pen), 
decagonal (m = 10, Dec), and dodecagonal (m = 12, Dod) local order is possible. We find 
that m-gons are not stable for other values of m. For n > 6 the vacancy is too big. The 
sequence of local orders as a function of tq is expected to look like in Fig. Wljo)- 

Next, we calculate the T = 0, p = phase diagram in the ro-e-cr^ parameter space 
via annealing simulations and structural relaxation. For details, we refer to [21], where 
the same procedure has been applied. The result is depicted in Fig. [21 A wide second 
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Figure 1: (a) The local lowest energy configurations in two-dimensional monodisperse 
systems are polygons with or without a central atom, (b) Expected phase behavior of 
two-scale potentials, tq is the position of the second minimum. 



Phase 


Dim 


Density 


Lattice constants 


Particles/u.c. 


Symmetry 


References 


Xi 


2 


0.76 


a = 4.24, a = 72° 


13 


cmm 


|24j 


Xi2 


2 


0.80 


a = 1.62, b = 3.08 


4 


Pgg 


m 


Sig 


2 


1.07 


a = 2.73 


8 


p4g 


miisiE 


Pen 


2 


0.90 


a = 2.62, b = 4.25 


10 


cmm 


mm 


Pen2 


2 


0.99 


a = 2.41 


5 


p31m 


m 


sig 


3 


1.01 


a = 3.86, b = 2.00 


30 


P42/mnm 


HI 



Table 1: The complex crystals observed in the LJG system. The nearest-neighbor distance 
is equal 1. 



minimum (a^ = 0.042) stabilizes mostly simple phases. There are two hexagonal lattices 
(Hex, Hex2), which are connected continuously along small e- values. A rapid increase of 
the lattice constant is observed at the dashed line. Additionally, the square lattice (Sqa) 
and a phase built from deformed pentagons and triangles (Pen2) are found. 

Further structures appear for cr^ = 0.02: a complicated phase with decagons (Xi), a 
related one with flattened honeycombs (Xi2), a phase with pentagons and hexagons (Pen), 
the two-dimensional sigma phase (Sig), and alternating rhombs (Rho). In the case of 
a"^ = 0.01, parallel rhombs (Rho2) are stable. We remark that the phase diagram differs 
slightly from the one published in [21] due to a higher precision of the relaxation: Sig2 and 
Sig3 are now unstable, whereas Xi2 and Pen2 appear as new phases. In total, nine ground 
states have been discovered, among them five complex crystals (Tab. [1]). 

The sequence of phases in Fig. |2] follows the local order analysis in Fig. [H Hence, the 
phase diagram can be understood as the result of stabilizing regular polygons. We assume 
that the behavior is universal in the sense that the phase diagram looks similar for most 
two-scale interaction potentials. 

It is interesting to study the ground state energy per particle, E, as a function of the 
potential parameters. In Fig. [3] the rescaled energy E'{rQ,e) = E/{e + 1) is plotted. The 




Figure 2: Phase diagram of the two-dimensional LJG system for three values of a^. For 
each crystal the decoration of a unit cell is shown. 
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Figure 3: Rescaled ground state energy E' for cr^ = 0.02. Maximum and minimum are 
located in the stability regions of Sqa. 

rescaling enforces convergence in the limit e ^ oo. For a given phase the energy has a 
minimum near tq defined by ([2]) and ([3]). Sometimes the phase is not stable at its minimum 
energy. As seen in Fig. [3l the stability regions of Pen/Pen2 and Xi/Xi2 are shifted to lower 
To values due to competition with Hex. 

An important aspect of the two-dimensional LJG system is the fact, that all phases can 
actually be grown as single crystals in simulations, if the cooling time is slow enough. When 
heating the system up, the complex crystal often reversibly transform into quasicrystals. 
This is the case for Xi [23], but also for Pen/Pen2 and Sig [25j, and means that these 
phases are not accessible from the melt, but have to form via a solid-solid transformation 
from the quasicrystal [26] . 

3 The sigma phase and a glass in three dimensions 

In this section, first results for the three-dimensional LJG system are presented. Its full 
T = 0, p = phase diagram is not known. Nevertheless as we show now, the LJG potential 
allows to grow at least one complex crystal in simulations. We ran 11 simulations with 
e = 1.8, (T^ = 0.02, and tq G [1.0,2.0], Aro = 0.1. Standard molecular dynamics (MD) 
with a Nose-Hoover thermo-/barostat and periodic boundary conditions were applied. The 
cubic simulation box contained 2744 particles with an initially amorphous configuration. 
We searched for temperatures close to, but below the melting point for rapid crystallization. 
We were successful to achieve crystallizations for all r^ values except 1.4. After a few 
10^ MD steps, six different crystals have been grown: fee (cF4, tq = 1.0), hep (hP2, 1.6), 





(b) 



(c) 



Figure 4: Single crystal of the sigma phase grown with molecular dynamics using 15 625 par- 
ticles, (a) Projection along the four-fold symmetry axis, (b) Diffraction image, (c) Fourier- 
filtered structure image with tiling superimposed. 



bcc (cI2, 1.1-1.2, 1.7-1.9), shl (hPl, 1.3), pzt (tP5, 2.0), sig (tP30, 1.5). Here, shl is the 
simple hexagonal lattice (prototype: 7-HgSn6_io, disordered), pzt the tetragonal perovskite 
structure (PbZrxTii-xOs), and sig the sigma-phase (/3-uranium) . 

It was observed that the particles cannot relax effectively with periodic boundaries, 
since point defects and dislocations do not heal out. Therefore, the simulations were 
repeated with open boundaries, i.e. the particles form a solid sphere floating in vacuum. 
The structures are the same, but this time single crystals were obtained. The possibility 
to simulate at zero pressure is an advantage of the LJG potential. For the Dzugutov 
potential, external pressure usually has to be applied, because of the repulsive bump in 
the interaction. 

For To = 1.5, the sigma phase was observed. To study its formation in detail, we 
simulated a large system of 15 625 particles at T = 1.7 (as usual ks = 1), which is about 
90% of the melting temperature. After less than 10^ MD steps, local crystallization started 
and a polycrystal with dodecagonal quasicrystallites appeared. Ordering proceeded in two 
steps: first to a bicrystal at 10^ steps and then to a single crystal at 2 • 10^ steps. At the 
same time, rearrangements within each crystallite transformed the quasicrystal state into 
the sigma phase. After 3-10^ steps the system was relaxed to T = 0. 

Fig. 111(a) shows the configuration of the particles projected along the four-fold axis. 
The width of the Bragg peaks in the diffraction image (Fig. IH^b)) is determined by the 
system size. By selecting only the twelve inner strong reflections for an inverse Fourier 
transform, we obtain a Fourier filtered structure image (Fig. IHc)). It allows to easily 
detect the underlying tiling. Cylindrical columns with pentagonal shape forming squares 
and triangles in a ratio 1:2 are characteristic for the sigma phase. 

We note that the sigma phase is unstable at T = 0, because fee has a lower energy than 
sig. To test the stability at T = 1.7, we initiated a simulation of a fee crystal in contact 



with a sigina crystal. Quickly, the whole system transformed into a single fee phase. Hence, 
the sigma phase is also not stabilized entropically. The reason for its appearance in our 
simulations is the huge nucleation radius of fee. The local order of sig is much closer to 
the melt. Further details will be reported elsewhere. 

In the case tq = 1.4, all attempts to crystallize the system failed. One reason is the 
comparably low melting temperature: the ground state fee melts around T = 0.9. In 
various MD runs over several 10^ steps in the range 0.5 < T < 0.9 no nucleation was 
observed. We conclude that the choice tq = 1.4 is interesting for studying a monatomic 
glass. Our results indicate that this LJG glass is more resistant against crystallization 
than the Dzugutov glass [27], which forms the dodecagonal quasicrystal rather quickly 

4 Discussion 

How does the LJG system compare to experiments? In metals, multi-body terms can only 
be neglected in a first approximation. Effective potentials often have Friedel oscillations, 
which are mimicked by double-wells. It should be kept in mind that fixed pair interactions 
between all atoms might be not applicable in complex phases, since the atoms are found 
in different local environments. Furthermore, the interaction is expected to change during 
crystallization. 

Isotropic pair potentials are more applicable to macromolecular self-assembly, because 
the molecules as a whole interact almost classically. It has been suggested [29J that their 
complex arrangements originates from the competition of two length scales, which appear 
due to soft repulsion and strong interpenetration. The micelles forming the dodecagonal 
quasicrystal in [15] have two natural length scales: the inner one corresponds to the back- 
bone of the dendrons and the outer one to the end of the tethered chains. Our simulations 
suggest that the ideal ratio of the scales is close to 2:3. 

The cylindrical phases (T-shaped molecules and ABC-star polymers) can be understood 
as two-dimensional tilings with an effective interaction within the plane. Although not 
isotropic anymore, the particles still have two length scales. The tilings observed so far 
[30} |3T] consist of hexagons, squares, and triangles only. If it will be possible to stabilize 
pentagons or decagons, then the phases Pen/Pen2, Xi/Xi2 might appear. 

We finish with a challenge: Can monodisperse icosahedral quasicrystals be grown in 
simulation or experiment? As of today, none have been found. 
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